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We present accurate numerical estimates for the correlation amplitudes of leading and main 
subleading terms of the two- and four-spin correlation functions in the one-dimensional spin-1/2 
XXZ model under a magnetic field. These data are obtained by fitting the correlation functions, 
computed numerically with the density-matrix renormalization-group method, to the corresponding 
correlation functions in the low-energy effective theory. For this purpose we have developed the 
Abelian bosonization approach to the spin chain under the open boundary conditions. We use 
the numerical data of the correlation amplitudes to quantitatively estimate spin gaps induced by a 
transverse staggered field and by exchange anisotropy. 
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I. INTRODUCTION 

The one-dimensional spin-1/2 XXZ model is a basic 
and well studied model in statistical physics. In some 
parameter range its ground state is critical and spin- 
spin correlation functions exhibit quasi-long-range order. 
The model represents a typical example of Tomonaga- 
Luttinger (TL) liquids in which elementary excitations 
are gapless collective modes rather than single-particle 
excitations, in contrast to three-dimensional systems. 
These features are uncovered by extensive theoretical 
studies over the years which have employed various pow- 
erful methods including the Bethe ansatz^ bosoniza- 
tion technique^ and the conformal field theory^ One- 
dimensional spin chains are also important from the ex- 
perimental viewpoint, as they are relevant to many quasi- 
one-dimensional magnets in which couplings along one 
direction are considerably stronger than those in the 
other two directions. 

In this paper, we discuss long-distance asymptotes of 
equal-time correlation functions in the spin-1/2 XXZ 
chain in a magnetic field. The Hamiltonian is 

Wo = J E ( W +1 + SfSf +1 + AS?S? +1 ) - H £ Sf , 



(1) 

where J > and Si — (Sf,Sf,S[-) is an S = 1/2 spin 
operator on the Zth site of the chain. The anisotropy 
parameter A is assumed to satisfy the inequality A > — 1 
such that the system is in the critical regime for a certain 
range of magnetic fields, i.e., < \H\ < H C 2 for —1 < 
A < 1 and H cl < \H\ < H c2 for A > 1 where H cl and 
H C 2 are the lower and upper critical fields, respectively. 
The spins partially polarized by the field H have a finite 
magnetization per site M, where —1/2 < M < 1/2. The 
low-energy excitations of TCq in the critical regime are 
free massless bosons. They are governed by the Gaussian 
theory 
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dx 
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where v is the renormalized spin-wave velocity and the 
bosonic fields 4>(x) and 4>{x) obey the commutation re- 
lation [<t>(x), 4>(y)] — — (i/2)[l + sgn(a; - y)]. We take 
the lattice spacing a = 1 and identify the site index I 
with the continuous space variable x. For a better de- 
scription of low-energy physics, one needs to add to 7io 
an irrelevant operator cos(20/i?): This operator becomes 
marginal at A — 1 and H = yielding logarithmic cor- 
rections, while it is relevant and opens a spin gap for 
A > 1 and < \H\ < H c \. We ignore this term in this 
paper, however. This may introduce systematic errors in 
our analysis for A > 1 and small M . The spin operators 
are related to the bosonic field a 3 i 4 i 5 as 
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where Q — 2itM is incommensurate wave number, and 
a n and b n are non-universal constants which depend on 
A and M. The TL liquid parameter (or the compacti- 
fication radius) R characterizes the asymptotic behavior 
of correlation functions and its exact value is readily ob- 
tained by solving the Bethe ansatz integral equationsALS 
Essentially all the low-energy properties of the spin chain 
JTJ follow from Eqs. ©-Q- For example, the ground- 
state two-spin correlation functions are found to decay 
algebraically^iiSi 
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The decay exponent ij — 2ttR 2 is a function of A and M. 
For -1 < A < 1, r) = 1 - cos _1 (A)/7r at M = 0, while 
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r] = 2 at M — » +0 for A > 1. As M increases, 77 varies 
monotonically and approaches the universal value 1/2 in 
the limit M — > 1/2. The amplitudes and are 
related to the coefficients a n and 6„. Recently Lukyanov 
and coworkerisiiii2iiii have obtained exact formulas of the 
correlation amplitudes in Eqs. (01 and © at M = and 
— 1 < A < 1. Nevertheless their values at M ^ are not 
known analytically. 

The aim of this paper is to numerically determine the 
non-universal coefficients a n and b n with high accuracy 
for arbitrary M. We first extend the bosonic representa- 
tion of the spin operators J3J and to a more general 
form of infinite series, and calculate the spin polarization 
(Sf) and the two-spin correlation functions (SfSf,) and 
(SfSf,) within the bosonization theory. The analytic for- 
mulas so obtained are used to fit numerical data which 
we compute by using the density-matrix renormalization 
group (DMRG) methodi^i 5 - The fitting parameters are 
the coefficients a n and b n . This scheme is basically the 
same as the one used in our previous studies jiLiSi and in 
this paper we provide more accurate numerical data of 
the coefficients a\, bo, and b\ for wider range of parame- 
ters. 

The information on the coefficients is crucial for quan- 
titative analysis on effects of perturbations to which 
the critical ground state of the spin-1/2 XXZ model 
has instabilities. Such perturbations include bond 
alternation(A£ii& next-nearest-neighbor coupling ji2i22i2i 
and transverse staggered magnetic field j^SiSi 2 ^ 5 . In the 
bosonization approach, a perturbation written in the 
original spin operators is translated into bosonic oper- 
ators through Eqs. 10 and Q. The impact of the per- 
turbation on the ground state, the vacuum state of the 
Gaussian model © , can be qualitatively estimated from 
scaling dimension of the perturbation operators, which 
depends only on the TL liquid parameter R. It is, how- 
ever, necessary to know the exact values of the non- 
universal coefficients a n and b n as well, to quantitatively 
analyze effects of the perturbation, e.g., to estimate mag- 
nitude of an energy gap induced by the perturbation or 
to compute correlation functions in the presence of the 
perturbation. As an example of such quantitative anal- 
ysis we calculate the spin gap induced by the staggered 
transverse field 22,23 - 24,25 from the numerical data of the 
coefficients. 

There are often cases when perturbations are compos- 
ite operators of two (or more) spins, such as exchange 
interactions. In this case one must consider fusion of two 
operators. In general, fusing two spin operators at short 
distances may generate operators which are not present 
in each single-spin operator which are fused. Conse- 
quently, the leading term of the correlator of the com- 
posite operator may be different from a product of the 
leading terms of the correlators of the single-spin oper- 
ators. We illustrate this by taking the nearest-neighbor 
coupling SfSf +1 as an example. We examine operators 
generated by the fusion of Sf and Sf, 1 . From a numerical 
fitting of four-spin correlation functions, we estimate the 



amplitude of a leading uniform term of the correlation 
function. The results are used to analyze the spin gap 
induced by the perturbation of exchange anisotropyiSSi 
The paper is organized as follows. We briefly review 
the Abelian bosonization approach to the model in 
the next section. The bosonic representation of spin op- 
erators in the form of infinite series is introduced in Sec. 

II A. We then use it to derive the analytic formulas of 
the two- and four-spin correlation functions in Sec. II B 
and C. In Sec. Ill, we present the numerical results on 
the correlation functions obtained from the DMRG cal- 
culation. The numerical data of amplitudes of the two- 
and four-spin correlation functions are presented in Sec. 

III A and B, respectively. We show that the numerical 
results are in excellent agreement with analytical predic- 
tions available for various limiting cases. Furthermore, 
we discuss the spin gaps induced by the perturbations 
of staggered transverse field and exchange anisotropy in 
Sec. Ill C. Finally, we summarize the results in Sec. IV. 
Appendix explains bosonization of spin operators. 



II. BOSONIZATION 

A. Spin operators 

In this section, we summarize the Abelian bosonization 
technique applied to the spin-1/2 XXZ chain in a mag- 
netic field. Let us first express the original spin operators 
Sf and Sf in terms of the bosonic fields. To this end, we 
follow and extend the scheme of Refs.y and|j. We begin 
with the bosonic representation of electron operators in 
the Hubbard model with the on-site repulsion U. At half 
filling we have a gapped charge mode and a gapless spin 
mode. After constructing spin operators from the elec- 
tron operators, we integrate out the gapped charge mode 
to obtain a bosonic representation of the spin operators 
in the Heisenberg chain. We then generalize the result 
to the XXZ case. A detailed derivation is presented in 
Appendix, and here we show only the final results, 
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where a n and b n are non-universal constants. They de- 
pend on the short-distance regularization of the Gaussian 
theory as well as on the parameters A and M. Equations 
© and are just the first few terms of Eqs. J7J and 
© with the shift of the field <p(x) — > 4>{x) + QRx. In 
principle the right hand side of Eqs. Q and JSJ| should 
contain descendant fields as well. It follows from Eq. © 
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Here we have used the commutator [4>(x), <f>(x)] = —i/2. 
An equivalent bosonic representation of the spin opera- 
tors is recently derived in Ref. from global symmetry 
analysis of the lattice and field operators. Our derivation 
in Appendix is complementary to Ref. Il3l 
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where (• • •) = (0| • • • |0) represents the expectation value 
in the lowest energy state in the sector with magnetiza- 
tion M. If we adopt the regularization 



oo -Afe 
dk-— 



(1 — cos kx) = lnx, 



(15) 



then the correlation amplitudes are related to the coeffi- 
cients in Eqs. Q and JHJ by 



B. Two-spin correlation functions 

The two-spin correlation functions are readily calcu- 
lated by the bosonization method. Let us first consider 
the thermodynamic limit where the correlation functions 
depend only on the distance between the two spins. In 
this case we expand the bosonic fields <f>(x) and <fi(x) as 
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where a>h and at are boson operators obeying the com- 
mutation relation [afc,^,] = 5{k — k'), A is a short- 
distance cutoff, and S is a positive infinitesimal. The 
fields 4>{x) and (j>(x) defined by Eqs. (|10|) and 111( 1 satisfy 
the commutation relation [4>(x), (j)(y)] = ~ i@(x — y) in 
the limit A — > 0, where Q(x) is the step function. The 
Gaussian Hamiltonian (J2J defined on the whole real x 
axis now reads 
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Thus, the ground state |0) of the Hamiltonian corre- 
sponds to the vacuum for the bosons a k , i.e., a k \0) = 0. 
Substituting Eqs. and JTTJ to Eqs. JTJ and ©, we 
obtain the equal-time two-point correlators in the ther- 
modynamic limit, 
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For M = and — 1 < A < 1, the values of the amplitudes 
Aq, Af, and A\ are obtained by analytica l 1 ^ 12 ' 1 ?! 2 ? and 
numerical^ methods. We note that in principle the long- 
distance expansion lfH5|l and iJTlj l should also include con- 
tributions from the descendants ignored in Eqs. J7J) and 
@ and those from the irrelevant operators discarded in 

As mentioned in Introduction, we determine the coef- 
ficients a n and b n by fitting numerically computed cor- 
relation functions to appropriate formulas. Since the 
numerical DMRG method works best for finite-size sys- 
tems with open boundaries, we need calculate the cor- 
relation functions under the open boundary conditions. 
For this purpose we employ the open-boundary bosoniza- 
tion scheme developed in Refs. and Suppose 
that the XX Z spin chain of our interest consists of L 
spins Si (I = 1,2, ... ,L). This is equivalent to assume 
So = Sl+i = 0. In the bosonic representation this 
amounts to have the Gaussian model @ defined in the 
finite region < x < L + 1 with Dirichlet boundary con- 
ditions at x = and x = L + 1. Our convention is that 



c/)(0) = 0, <j>(L + 1) = 2-kRLM. 



(17) 



These boundary conditions are consistent with Eqs. {jj- 
(jnj. Once we fix the boundary conditions (jT7|l and the 
regularization scheme, the coefficients a n are uniquely 
determined whereas the coefficients b n are determined 
only up to a phase factor. Instead of Eqs. ljTU|) and ltTTJl . 
we now have the mode expansion&ii2ii& 
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where q n — mi/(L + 1), [(j>o, 4>o] = h and a n and Sj, are 
boson operators satisfying [a^o;^,] = The Gaus- 
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sian model becomes 
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2{L + 1) 24{L 
The lowest-energy state of the spin chain Q with mag- 



netization M corresponds to a vacuum of bosons a n |0) = 
with 0o 1 0) = 27ri?LM|0). To calculate the zero- 
temperature correlation functions, we substitute Eqs. 
lfT5)l and l(T$|) into Eqs. J7J) and © and take average with 
respect to the state |0). (The readers interested in the 
detailed calculation should refer to Ref. 0.) We obtain 
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and the sum i s taken over odd n only. We have used 
the regularization J2 n =i\^~ cos ( < l'n x )}/ n = m [/i( x )] as m 
our previous studies^S^ The factor s(n,n' '',1,1') in Eq. 
(EH is 

s(n, n'; Z, Z') 

' (-l)(™+ 1 )'+("' + 1 ) i ' + (™+«')/2 
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In the thermodynamic limit {L — > oo with |Z — L/2| <C 
and \V — L/2\ <C L) the correlators and reduce 
to Eqs. lO and 1(3. 



C. Four-spin correlation functions 

In this subsection, we discuss fusion of two operators 
taking Sf Sf +1 as an example. To find bosonic represen- 



tation of the composite operator SfSf +1 , we need oper- 
ator product expansion of the operators in Eq. 10 . We 
explicitly write down the product of Sf and Sf +1 as 



s T s i+i = Jq E E (26) 

where 

X({e i( e-,ni};Z) = b ni b n2 t{{ei, e-, n,; Z}) 

x exp |i27ri?[ei0(Z) + e 2 0(Z + 1)]} 

x exp |-^Kei0(Z) + n 2 e' 2 (j){l + 1)] 

(27) 

and t{{ci, n^; Z}) is defined by i = — 1 {ni,n,2 = even), 
t = -eie 2 e' 1 £2 (n l7 n 2 = odd), t = — i{— l)'e 2 e' 2 (rii = 
even and n 2 = odd), t — i{—lye\e'i (ni = odd and n 2 = 
even). To find the first few leading operators in the ex- 
pansion (|26|l . we make each exponential operator in Eq. 
(|27|l in normal order and expand the fields </>(Z + 1) and 
4>{l + 1) as <j){l + 1) = <f>{l) + d(f>{l)/dl + ■•-. It turns 
out that the leading operators in Eq. 112611 come from the 
following three contributions: 

(i) ei + e 2 = and nxe{ + n2€ 2 — 0. In this case X 
is expanded as a sum of a constant (identity operator) 
term, d<j>/dx and d(p/dx with scaling dimension 1, and 
higher-order terms. 

(ii) t\ + e 2 = and nie^ + n 2 e 2 = ±1. The leading 
operators in the expansion of X are (—1)' exp[±i(p{l) / R], 
whose scaling dimension is 1/ {2rf) . 
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(iii) e± + t 2 = ±2 and nie[ + n 2 e' 2 = 0- The leading 
operators in this case are exp[±i47ri?0(/)], which have 
dimension 2rj. 

We thus find that the operator SfSf +1 has the eXpan- 
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Apart from the constant term Co, the leading term in 
SfSf+l is the oscillating one with the coefficients c\ and 
c[ if A > (77 > 1/2), while it is the non-oscillating term 
with Cq and c' Q ' if A < (7/ < 1/2). We note that, for both 
signs of A, the leading term in the expansion of SfSf +1 is 
not simply given by the product of the leading operators 
in Sf and Sf +1 . In fact, the higher-order terms with large 
n in Eq. JJ]J, which can be ignored in the calculation of 
long-distance asymptotes of two-point correlation func- 
tions, give contributions to the leading operator in the 
operator product expansion of SfSf +1 . The same ob- 
servation can be made for the other subleading terms in 
SfSf+i' This result illustrates that, when one fuses two 
(or more) operators at short distances, information on 
only a few leading terms of each operator is, in general, 
not enough to determine the leading terms in the bosonic 
representation of the fused operator. Instead, one must 
survey contributions from all higher order terms. 

Calculating the coefficients c„ from the coefficients b n 
is hardly possible not only because each c„ has contri- 
butions from infinitely many 6„s but also because of the 
ambiguity due to the short-distance cutoff. It is thus 
more practical to estimate the c„s from the four-spin cor- 
relation function (Sf Sf +1 Sf, Sf, +1 ) and its variants. From 
Eq. H28fl we know the asymptotic behavior of the four- 
spin correlation functions, 
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where : O := O — (O). The correlation amplitudes are 



related to the coefficients 



by B = 4, B 1 = (c{ 
(cq 2 + < 2 )/2^, and B 2 = c|/2. 
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FIG. 1: (a) (-ljI'-'^Sf SJ5> versus r = \l - l'\, (b) \(Sf Sf,)\ 
versus r, (c) KSf)! versus I for A = 0.5, 0, and —0.5 and 
M = 0.25. The open symbols are the DMRG data while the 
small dots are the results of fitting. 
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III. NUMERICAL RESULTS 

In this section, we present numerical results on the cor- 
relation amplitudes in the two- and four-spin correlation 
functions obtained from the DMRG calculationiiii^ We 
calculate the spin polarization (Sf), the two-spin corre- 
lation functions (SfSf,) and (SfSf,), and the four-spin 
correlation functions (:(S~ h Sf +1 + SfSfl-^): :(SyS~7 +1 + 

S F S t/+i) : > and ( S t S l~+i S l/ S l/+i) in the °P en XXZ 
chain of L — 200 sites. The correlation functions are 

calculated for I = r — r/2 and I' = r + r/2, where r 
represent the center position of the chain, r — L/2 (for 
even r) or r = [L + 1) /2 (for odd r). The numerical cal- 
culation is done using the finite system algorithm, and 
the number of kept states m is up to 200. We estimate 
the numerical error due to the DMRG truncation from 
difference between the data computed with m = 200 and 
those with m = 150. The estimated errors for the spin 
polarization, two- and four-spin correlation functions are 
typically less than 10 -7 , 10~ 6 , and 10~ 6 , respectively, 
and sufficiently small for accurate estimation of the am- 
plitudes. 

A. Amplitudes of two-spin correlation functions 

First, we show the results on the spin polarization 
(Sf) and the two-spin correlation functions (SfSf,) and 
(SfSf,). Since the nth order terms in Eqs. Q and © 
contribute to the correlators less and less for large n, we 
may neglect the higher order terms with n > 2 in the 
fitting procedure. That is to say, we fit the DMRG data 
to the analytic formulas l|21[l -H23 [) setting a n = b n = 
for n > 2 and taking bo, b\, and a\ as fitting parame- 
ters. We note that this scheme for determining the coef- 
ficients is basically the same as those used in our previous 
studieSfi2ii& in oneiS of which the numerical data of Af 
and Aq are reported for several typical values of M and 
< A < 1. However, in that workiS the decay exponent 
•q as well as the coefficients was taken as a fitting param- 
eter, and this could cause small but avoidable errors in 



the estimates of the coefficients. In the present work, we 
use the exact value of rj obtained from the Bethe ansatz 
solutions. We therefore believe that the estimates of the 
coefficients presented here are even more accurate than 
the previous ones. 

Figure shows the two-spin correlation functions 
(SfSf,), (SfSf), and the spin polarization (Sf) at A = 
0.5, 0, -0.5 and M = 0.25. The DMRG data and the 
fitting results are plotted with the open and solid sym- 
bols, respectively. The excellent agreement between 
them demonstrates that the fitting procedure works ex- 
tremely well at least for the parameters used in Fig. Q 
To determine the correlation amplitudes, we perform the 
fitting for the data of several ranges, 20 < r < 140, 
20 < r < 180, 60 < r < 140, and 60 < r < 180 for 
the two-spin correlation functions, and 20 < I < 180, 
40 < I < 160, and 60 < I < 140 for the spin polar- 
ization. We then take the mean and the variance of the 
fitting results as the estimate and the error of the numer- 
ical values of the amplitudes, respectively. The results 
of = bl/2, Af = 6?/4 estimated from (SfSf,) and 
A\ = af/2 estimated from (Sf) are plotted in Fig. [5] as 
functions of M for several typical values of A. As shown 
in the figure, the amplitudes take non-universal values at 
M = Q and vary smoothly as M increases. We have con- 
firmed that the data of Af estimated from the correlation 
function (SfSf) coincide with those obtained from (Sf) 
within error bars. We note that, for small A < —0.8, 
the accuracy of the estimated amplitudes Af and Af be- 
comes considerably poor. This difficulty might be due to 
the fact that, for this parameter range, the subleading 
terms cx Af,Af of the correlation functions become con- 
siderably smaller than the non-oscillating leading terms. 
Further studies with a more elaborated scheme will be 
required for accurate estimation of Af' z in this case. 

To further examine the accuracy of the numerical es- 
timates obtained above, we compare them with exact 
results which are available for some limiting cases. 

At M = 0, the correlation amplitudes Aq, Af and A\ 
for —1 < A < 1 are analytically calculated>ii*i2ii£ 



An = 



Af 



A 7 ; 



8(1-77)2 

1 

277(1 - 77) 



re n 1 

V^r(^) 
r( a \ 

2^r( 5(I ^ y ) J 



exp 



dt 1 sinh(77i) 



t \sinh(t) cosh[(l — rj)t] 



(32) 



TV U. \ 



(!->?)' 



exp 



exp 



dt 



dt 



cosh(2ni)e- 2t - 1 



o t \2smh(r)t) sinh(t) cosh[(l — rf)t] sinh^i) 77 



" + V 2 * 



sinh[(2r? - l)i] 2n - 1 



-21 



t \sinh(7yt) cosh[(l — rj)t] rj 
I 



(33) 
(34) 



where T(x) is the Gamma function. These equations were previously confirmed by numerical calculations^*^^ 
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FIG. 2: Amplitudes of the two-spin correlation functions as functions of M for several typical values of A; (a) Aq, (b) Af, and 
(c) Af. Aq and Af shown in (a) and (b) are estimated from (SfSfr). A\ for M > in (c) are estimated from (Sf) while A\ at 
M = are from (S* Sj/). The solid curves in (a) and (b) represent Eqs. 135H and 1361 . respectively, while the solid line in (c) 
represents the exact result Af = 1/(2tt 2 ) valid at A = 0. 



and here we have found that the numerical estimates of 
the present work with higher accuracy are even in better 
agreement with the above exact formulas for — 0.8 < A < 
0.8. 

In the saturation limit M — > 1/2, exact asymptotic 
form of the correlation (SfSf,) can be obtained from 
known exact results on the hard-core boson modeli22iiffii2i 
It follows that near the saturation limit the amplitudes 
Aq and Af should behave as 

A ° = 2^ (J"")'' (35) 

^•-is?w(i -")"'■ <36) 

where p is a universal constant related to Glaisher's con- 
stant A by p = ne 1/2 2-^ 3 A- e = 0.92418 • • •. We can see 
in Fig. |21 (a) and (b) that, for all As shown, the numer- 
ical data of Aq and Af approach the predicted behavior 
shown by the solid curves. Note that there is no free 
parameter in the theoretical predictions 1)35(1 and l|Ho[l. 

As for the {SfSf,) correlation, it is expected from Eq. 
()14|l that the amplitude Af should converge to a universal 
value l/(27r 2 ) at M — * 1/2 for arbitrary A, since in this 
limit r) becomes 1/2 and the correlator must take the 
constant value 1/4. Furthermore, for A = one can 
easily calculate (Sf Sf,} exactly using the Jordan- Wigner 
transformation to find Af = l/(27r 2 ) for arbitrary M. 
We clearly see that the numerical data of Af in Fig. 0(c) 
agree with these predictions. 

From these observations, we conclude that our esti- 
mates for the correlation amplitudes Aq, Af, and Af 



are pretty accurate, except for the following parameter 
regime: (i) A < -0.8 and (ii) A > 0.8 and M = 0. In the 
latter regime the leading irrelevant operator cos(20/i?) 
is no longer negligible. We point out that diverging be- 
havior at M — 0, due to the presence of the (almost) 
marginal operator, can be clearly seen in the data of Aq 
and Af for A > 0.8. 

We summarize the numerical data of bo, b\, and a\ in 
Table [I] It is important to note that the sign of a\ is 
uniquely determined once the boundary conditions l|17|) 
are fixed. On the other hand, bo and b\ have the ambi- 
guity of a common phase factor. In the Table, we show 
the data with fixing bo to be real and positive. 

B. Amplitudes of four-spin correlation functions 

Next we discuss numerical results of the four-spin cor- 
relation functions. Figure [3] shows the numerical data of 
the four-spin correlation functions ( : (S^ S^ +1 + S' i ~<S , ; ~|j_ 1 ) : 
: + S-S+ +1 ) : ) and (S+S+ for A = 

0.5, 0, and —0.5 and M = 0.05. We see that the corre- 
lator (:(S+Sr +l + SfS+J. :(S+S^ +1 + S^S+ +1 ):) shows 
power-law decay with the exponent equal to either 1/rj 
(for A > 0) or 2 (for A < 0), while (SfSf +1 S^Sy +1 ) 
always decays with the exponent 4?/. We have found the 
same behavior for other values of A and M. These results 
are consistent with Eqs. H29fl -131 [) . 

To estimate the amplitude B2 — c|/2, we fit the nu- 
merical data of (S'f S~ l * +1 Sf SjJ +1 ) to its analytic formula 
for finite chains with open boundaries, 

(^l S l+1 S { , 5)/+i) 
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A=0.5, M=0.05, L=200 




10 r 

A=0.0, M=0.05, L=200 




O : |<:(S S~+S"S ): :(S S~+S~S ):>| 
□ ■ : <S + S + S~S~> 



10 10 



A=-0.5, M=0.05, L=200 




FIG. 3: Four-spin correlation functions |{: (S^S l+1 + 
S~S+ +1 ) :: (S+S~ +1 + S~S+ +1 ) :)| (open circles) and 
(5 ; + S^jSjT S^ +1 ) (open squares) versus r = \l — l'\ for M = 
0.05 and (a) A = 0.5, (b) A = 0, and (c) A = -0.5. The 
solid, dotted, and dashed lines correspond, respectively, to 
the algebraic decay of r~ 4v , r _1// '', and r~ 2 . The fitting re- 
sults for {Sf Sf +1 S^ S^ +1 ) using Eq. 1371 are plotted by solid 
squares. 



B 2 

0.015 



0.01 



0.005 




FIG. 4: Amplitude B2 of the four-spin correlation func- 
tion {SfSi +x SyS^ +1 ) as functions of M for several typ- 
ical values of A. The solid curve represent the relation 
Bi = [1 + cos(27rA/)]/167r 2 expected for A = 0. Inset : Nu- 
merical estimates of the amplitude B2 at M — (open circles) 
and the analytical prediction Eq. Q380 (solid curve). 



(2c 2 ) 2 (exp iAnR4>(l) exp -UnR4>(l') ) 



8B 2 



/ 2 „(2Z)/ 2 „(2H 

hvii + iOMi-i'V 



(37) 



where we used Eq. I|19|) ■ We see in Fig. [3] that the data 
of (S^Sj, 1 SfS^7 +1 ) are fitted by the formula extremely 
well. The estimated values of B2 from the fitting proce- 
dure are shown in Fig.0] We see that for each A the am- 
plitude takes a non-universal value at M = 0, decreases 
monotonically as M increases, and vanishes eventually 
at M -> 1/2. We note that for A = the analytic 
form of (Sf S^Sy S^ +1 ) can be easily calculated, yield- 
ing B 2 = [1 + cos(27rM)]/167r 2 . Figure H shows that the 
numerical data for A = are in good agreement with the 
formula. We also compare the numerical estimate of B2 
at M = and — 1 < A < 1 with the exact formula of B2 
derived recently by Lukyanov and Terras^ 



B 2 = 



2 3+4r, 7r 2+2r )( - 1 _ 77 )2 



r( 



- 4— 4ri 



L \ 2-2-q) 



(38) 



Here again we find good agreement between the exact 
result and our numerical data. This is another evidence 
that our estimates are highly reliable. At M = and A > 
1, B 2 shows a diverging behavior due to the marginal 
operator, suggesting the break down of our analysis. 

Unfortunately, we cannot achieve a precise estimation 
of the amplitude B\ of the leading oscillating term in 
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(: (S+Si +1 + S^S+ +1 ) : : (S+5,7 +1 + SfS+ fl ) :) due to the 
presence of subleading terms which give sizable contribu- 
tions to the correlation function. This issue of estimating 
B\ is left for future studies. 



C. Spin gap 

The data of the correlation amplitudes obtained in the 
preceding subsections is useful for analyzing effects of 
perturbations of single-spin type and exchange-coupling 
type in the bosonization framework. To illustrate how 
this scheme works, we compute spin gaps induced by such 
perturbations to the Hamiltonian |JT]|. 

As an example of the perturbation of the single-spin 
type, we consider effects of the staggered transverse field. 
The perturbation to the Hamiltonian JTJ is given by 



H' = -h s J2(-l) l Sf 



(39) 



It has been shown that Ti' induces a spin gap. This 
field-induced gap is believed to be the origin of the spin- 
gap behavior observed in Cu benzoate^SiSiSiiSSi^ and 
Yb4As3^iiMi2£ under a uniform field, in which the stag- 
gered field emerges due to the alternating g-tensor and 
the Dzyaloshinskii-Moriya interaction. In these materi- 
als, exchange anisotropy is negligibly small and the stag- 
gered transverse field h s is proportional to the uniform 
field H . Thus, we may set A = 1 and h s = jH, where 7 
is a constant specific to each material. 

In the bosonization scheme, the leading uniform term 
of the perturbing Hamiltonian Ti', which is responsible 
for opening the gap, is a cosine term, 



Ti' = -h s b {H) / cos 2irR4>(x) 



dx. 



(40) 



Effects of the perturbation have been studied in the 
literatur o 22 i 23 i 24 i 25 i 36 i 37 and the induced spin gap is given 



E g 2v(H) T(g4 



2, t ' 



J 



nb (H) T(- 



2v(H) rm j 



(41) 



Hence, the field-dependence of the spin gap is evaluated 
by substituting our numerical estimates for bo(H) as well 
as the exact values of v(H) and rj(H) into Eq. (|41ll . The 
result is shown in Fig. 0for several typical values of 7. 
It reproduces the peculiar independence of the spin gap 
observed in experiments, E g ~ H 2 / 3 for small H. 

Next, we consider the spin gap induced by the pertur- 
bation of exchange anisotropy, 



n" = 



a) sfSu 



1 



(42) 



to the Heisenberg chain Tio(A = 1), where 1 — A <C 1. 
Note that, by rotating the system around the y-axis, the 



E n /J 
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:y=0.1 





: ^0.05 


• 


: y=0.01 
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/ ..-a 
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-C 



0.5 



H/J 



FIG. 5: Field-dependence of the spin gap in the Heisenberg 
chain H,o(A = 1) in a uniform field H and a staggered trans- 
verse field h s — ■yH . The dotted curves represent the expected 
power-law behavior of the gap, E g oc H 2 ^ 3 . 



whole Hamiltonian is rewritten as 

Ho + n" = jJ2 (sfsf +1 + SfSf +1 + AS?S[ +1 ] 

-H^Sf. (43) 
1 

Hence, the system can be also viewed as the 5=1/2 
XXZ chain with anisotropy A in a uniform transverse 
field H . As we have seen in Eq. Q28[l . apart from a con- 
stant, the leading operators in the bosonized Ti." are 

and -jj. Since their main effect is just a small renor- 
malization of the TL liquid parameter i?, we can neglect 
these operators in lowest order in 1 — A. We thus find 
that the dominant component which is responsible for 
opening the gap is the cosine term with coefficient C2, 



U" = -(1 - A)c 2 (H) / cos 4nR4>{x) dx 



(44) 



The effect of this perturbation has been studied^SiS 7 . and 
the spin gap for H ^> (1 — A) is found to be2i 



E g _ 2v(H) r(g^) 



J 



r( 



tt(1- A)c 2 (H)T(l-r)) 



2v(H) 



(45) 

We show in Fig. El the field dependence of the spin gap 
for several typical values of A. As shown in the figure, 
the spin gap opens very slowly with H and closes at the 
saturation field H = 2 J, reflecting the fact that the co- 
efficient C2 vanishes as M — > 1/2. In contrast to the case 
of the staggered transverse field, the gap induced by the 
exchange anisotropy is extremely small. This result is 
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0.001 



FIG. 6: Field-dependence of the spin gap in the Heisen- 
berg chain 7io(A — 1) with the perturbation of transverse 
exchange anisotropy (1 — A). 



consistent with the observation of the recent numerical 
study2& which finds no substantial gap; the spin gap is 
too small to be detected by the numerical study on finite- 
size systems. We note that the bosonization scheme with 
our estimates of ci is, at present, the only way to get 
reliable quantitative results on the spin gap behavior. 



IV. SUMMARY 

In this paper, we have studied the ground-state corre- 
lation functions in the 5=1/2 XX Z chain in a mag- 
netic field. With the bosonic representation of the spin 
operators Sf and Sf for open boundary conditions, we 
have calculated the two- and four-spin correlation func- 
tions analytically within the effective low-energy theory. 
We have also calculated the correlation functions numer- 
ically using the DMRG method, and estimated correla- 
tion amplitudes of the first few leading terms by fitting 
the numerical results to the analytic formulas. We have 
thus obtained precise data of non-universal coefficients 
appearing in the bosonic representation of lattice spin 
operators. Excellent agreement is found in the compar- 
ison of the numerical data with the exact known results 
in various limiting cases. 

We believe that the data of the correlation amplitudes 
presented in this work is suitable for quantitatively study- 
ing low-energy properties of perturbed spin chains within 
the bosonization method. Indeed, it has been shown in 
Ref. |3^ that the data of a\ and b can be used success- 
fully to explain quantitative features of the dynamical 
spin structure factor in Cu benzoate. We hope that the 
data will be applied to a wider variety of problems in 
one-dimensional spin systems. 
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APPENDIX A: BOSONIC REPRESENTATION 
OF SPIN OPERATORS 

In this appendix, we briefly overview the derivation of 
the bosonic representation of the spin operators J7J and 
(JSJ in the spin-1/2 XX Z chain. We basically follow the 
strategy of Refs. y and0; We first bosonize the repulsive 
Hubbard chain at half filling, in which charge excitations 
have a Mott-Hubbard gap. We then obtain the spin op- 
erators by throwing the gapped charge mode away, and 
generalize the results to the anisotropic XX Z case. 

Let us begin with the bosonic representation of electron 
operators (a =|, |), 




(Al) 
(A2) 

(A3) 



where Eire Klein factors obeying /V'} = 2<5 CT Cr ' 
and the bosonic fields ip obey the commutation relations, 

[<PR, a (x), <pR, a > (y)\ = (i/4)<Wsgn(x - y), (A4) 
kdWi^ff'fe)] = -(i/^Sa.a'Sgnix - y), (A5) 
k.W.Pvfe')] = -(i/4)S*,*>- (A6) 

The Fermi wavenumbers kpo are functions of the mag- 
netization M, k n = n(l/2 + M) and k Fl = n(l/2-M). 
We introduce fields <fi and <j) given by 

(j><j(x) = tp Lt<7 (x) + lp R<(J (x), (A7) 
4> a (x) = fL,a(x) ~ tpR,<r(x), (A8) 

which satisfy [<j> a (x), <j> a {y)] = -(i/2)[l + sgn(x-y)]5 CT)CT /. 
It then follows that the electron density becomes 

= i>R. a i>R,<T + i>L,^L,a 



IT 



1 



\fii dx 

The uniform charge and spin densities are 

1 dch 



(pt + pi)= 2 + ; 



^TT dx 

1 d(f> s 
/2tt dx 



(A9) 

(A10) 
(All) 
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where the charge and spin fields are defined by 



1 

71 
i 

71 



= - 7 =(0 T -0 i ).(Al3) 



The charge mode in the Hubbard chain is gapped 
at half filling when the on-site interaction is repulsive, 
U > 0. The charge gap is generated by the Umklapp 
scattering term, 



U 



Vr,\ ( x ) ^l.T ( x ) Vr,i ( x )^l,i (x) 

2U 
~ ~ (27ra) ; 



■ cos(v87T0, 



(A14) 



which pins the charge field at <f> c — n^ir/2 (n: integer). 
At low energies we may treat the field as a classical num- 
ber, i.e., 



(A15) 



cos(v87r0 c ) = C, 
sin(v / 87r0c) = 0, 



where C is a positive non- universal constant. At this 
point, following Ref. Isil we modify Eqs. (|A2|) and (|A3|I 
to 



'2na 

K. 



n=0 

oo 



'2na 



- i ( 2 n + 1 ) ( fc F a x + y^r 4> a ) — i V^4><r 



Using the equations above, one can derive the bosonic 
representation for the spin operators. The z-component 
of the spin operator is given by 



S'(x) 
_ 1 

~~ 2 
_ 1 

~ 2 



- - 1pl d (x)lpR,l(x) 



M + 



1 d<fi s 



2n dx 

oc 

E a 2n +i(-l) x sin (2ti+1)(2ttMx + V^k4> s ) , 

(A16) 



n=0 



where d2n+i is a non-universal constant. Here we must 
recall that Eq. i|A16[) is obtained from the Hubbard chain 
at half filling, whose low-energy effective spin Hamilto- 
nian is nothing but the antiferromagnetic Heisenberg spin 
chain, in which R = l/>/2n. To generalize the result to 
the XX Z chain, what one needs to do is replacing ^/27r<j) s 
with 4> S /R- We also define <fi(x) = 4> s (x) + 2ttRMx to ob- 
tain 



S\x) 



1 d(j>(x) 
2-kR dx 

OO 

-y^Q2»+i(-l) a: sin 



(2n+l) 



4>{x) 



R 



(A17) 



Similarly, the operator S + in the antiferromagnetic 
Heisenberg chain is given by 

S + (x) 

oo 

= e iV ^*° {b2n(-l) x cos 2u(2ttMx + V2^> s ) 

+ b 2n +i sin \(2n + 1){2ttMx + Vzirfa) } . 

(A18) 

By generalizing the equation to the XX Z case, we arrive 
at the final formula, 



S + (x) 



n=0 



2n- 



(x) 



R 



b 2n +i sin 



(2n + l) 



4> s (x) 



R 



(A19) 



where we have replaced \j2 , k^> s with 2irR(f> 
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TABLE I: Amplitudes (a) bo, and (b) 61, and (c) a± as functions of the magnetization M for several typical values of A. The 
values of a\ in (c) for M > are estimated from the fitting of (Sf) while those for M — are from (Sf Sf/). The figures in 
parenthesis indicate the estimated error on the last quoted digits. The error of the data without parenthesis is estimated to be 
less than 10" 4 . The exact values of Eqs. fgJ-JJD for M = are also listed. 
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0.5170 


0.5210 


4949H "1 


0.4984 


5096 


5068 


0.5109 


0.5149 


0.5188 




0.4968 


5008 


0.5048 


5088 


0.5127 


0.5165 


491 3f1 "l 


0.4951 


0.4989 


0.5027 


5065 


0.5103 


0.5139 


4896H ^ 


4Q39H 1 


0.4968 


5005 


0.5041 


0.5077 


0.5113 


4877(4 "1 

L/. 4:0 1 1 III 


0.4911(1) 


0.4946 


0.4981 


0.5016 


5050 


0.5084 


4857(1 "1 

\J • T:<_> !j 1 111 


0.4890 


0.4922 


0.4955 


0.4988 


0.5021 


0.5054 


4836H ^ 


4866H 1 


0.4897 


0.4928 


0.4959 


0.4991 


5099 


481 9(4 "1 

\J • 4:0 _1_ Ju \ -L ) 


4840H 


0.4869 


0.4899 


0.4929 


0.4958 


0.4987 


4787(4 "1 


481 3f1 1 


4840 


4868 


0.4896 


0.4924 


0.4951 


U.4: 1 OiJ 1 -L ^ 


fl 4784('1 


fl 480Q 
u.iouy 




fl 48R1 


fl 4887 


n 4Q1 3 


U.4: 1 OU^ -L } 




D 477fi 




0.4824 


fl 484Q 


fl 487^ 


U. 4:U£70 1 -L f 


0.4719(1) 


0.4740 


0.4762 


0.4785 


4808 

yj . 4:0*J0 


0.4830 


4663 fl") 


4683 fl 1 ) 


0.4702 


0.4723 


0.4743 


0.4764 


0.4785 




n 4(S44('1 "1 






fl 4(SQQ 


fl 471 8 


fl 47^8 




n 4fin9('i 


fl 4R1 R( 1 "1 




fl 4(S c i9 


fl 4R7D 


fl 4R87 


4543f1 1 


4557H 1 


0.4572 


0.4587 


0.4602 


0.4618 


0.4634 


44Qfif1 ^ 




0.4522 




fl 4 c i4Q 


fl d^P,^ 


n 4^77 


U.4:4:4:Ul _1_ 1 




fl 44(SS 


fl 44RD 


0.4492 


fi 4 c in c i 


n 4^1 7 


4^Q1 fl ^ 
u.ioyi^ J. ^ 


4401 f 1 ^ 


0.4411 


0.4421 


0.4432 


0.4443 


n dd^ 


U.4:OOZr^ -L ^ 


4^41 f 1 ^ 

U. 4:04: _L I _L 1 


n 4S4Q 


\J .1000 


fl 4^(S7 


fl 4^7(S 


n 4^8^ 


U.4:ZrUO^ ± ^ 




fl 49£9f 1 "1 


fl 498Q 


fl 49Q7 


fl 4^0^ 


n 4^1 ^ 


41QSf1 ^ 


42D4('1 1 


4210('1 , 1 


0.4216 


0.4222 


0.4228 


0.4235 


0.4122(1) 


0.4127 


0.4131 


0.4136 


0.4141 


0.4146 


0.4151 


403QH "1 


0.4042 


4046 

(J . 4:U4U 


0.4049 


0.4053 


0.4056 


4060 

yj .'rtuyjyj 


0.3947(1) 


0.3949(1) 


0.3952 


0.3954 


0.3957 


0.3959 


0.3961 


0.3845 


0.3847 


0.3848 


0.3850 


0.3851 


0.3853 


0.3854 


0.3731 


0.3732 


0.3733 


0.3733 


0.3734 


0.3734 


0.3735 


0.3602 


0.3602 


0.3602 


0.3602 


0.3602 


0.3602 


0.3601 


0.3455 


0.3454 


0.3454 


0.3453 


0.3453 


0.3452 


0.3450 


0.3283 


0.3282 


0.3281 


0.3280 


0.3279 


0.3278 


0.3276 


0.3072(3) 


0.3071(3) 


0.3070(3) 


0.3069(3) 


0.3068(2) 


0.3066(2) 


0.3064(2) 


0.2800(1) 


0.2800(1) 


0.2799(1) 


0.2798(1) 


0.2796(1) 


0.2795(1) 


0.2793(2) 


0.241(4) 


0.241(4) 


0.241(4) 


0.241(4) 


0.240(4) 


0.240(4) 


0.240(4) 



O.OO(exact) 



0.5141 



0.5079 



0.00 
0.01 
0.02 
0.03 
0.04 
0.05 
0.06 
0.07 
0.08 
0.09 
0.10 
0.11 
0.12 
0.13 
0.14 
0.15 
0.16 
0.17 
0.18 
0.19 
0.20 
0.21 
0.22 
0.23 
0.24 
0.25 
0.26 
0.27 
0.28 
0.29 
0.30 
0.31 
0.32 
0.33 
0.34 
0.35 
0.36 
0.37 
0.38 
0.39 
0.40 
0.41 
0.42 
0.43 
0.44 
0.45 
0.46 
0.47 
0.48 
0.49 



0.5006 
0.4937 
0.4904 
0.4883 
0.4871 
0.4863 
0.4857 
0.4853 
0.4850 
0.4847 
0.4844 
0.4842 
0.4839 
0.4835 
0.4830 
0.4825 
0.4818 
0.4810 
0.4801 
0.4790 
0.4778 
0.4764 
0.4749 
0.4731 
0.4712 
0.4690 
0.4666 
0.4639 
0.4610 
0.4578 
0.4543 
0.4504 
0.4462 
0.4416 
0.4365 
0.4310 
0.4250 
0.4183 
0.4110 
0.4029 
0.3940 
0.3841 

0.3728 
0.3601 

0.3455 
0.3284 
0.3074 
0.2802 
0.241(4) 



0.521(4) 
0.5022 
0.4983 
0.4960 



0.4945 

0.4934 

0.4926 

0.4919 

0.4914 

0.4909 

0.4904 

0.4899 

0.4895 

0.4889 

0.4883 

0.4876 

0.4868 

0.4859 

0.4849 

0.4838 

0.4825 

0.4810 

0.4795 

0.4777 

0.4757 

0.4736 

0.4713 

0.4686 

0.4658 

0.4627 

0.4593 

0.4556 

0.4517 

0.4473 

0.4425 

0.4374 

0.4317 

0.4255 

0.4188 

0.4114 

0.4032 

0.3942 

0.3842 
0.3729 
0.3601 
0.3455 

0.3284(1) 

0.3073(3) 

0.2802(1) 

0.241(4) 



0.512(2) 
0.5050(1) 
0.5029 
0.5014 
0.5003 
0.4994 
0.4986 
0.4979 
0.4973 
0.4967 



0.4960 

0.4954 

0.4947 

0.4939 

0.4931 

0.4922 

0.4912 

0.4901 

0.4888 

0.4875 

0.4860 

0.4843 

0.4826 

0.4806 

0.4784 

0.4761 

0.4736 

0.4708 

0.4678 

0.4645 

0.4609 

0.4571 

0.4529 

0.4484 

0.4435 

0.4382 

0.4325 

0.4262 

0.4193 

0.4118 

0.4035 

0.3944 
0.3843 
0.3730 
0.3602 
0.3455 

0.3283(1) 

0.3073(3) 

0.2801(1) 

0.241(4) 
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(a) b (0 > A > -0.6, A = 1.5, 2.0) 



M 


A = 0.0 


A = -0.1 


A = -0.2 


A = -0.3 


A = -0.4 


A = -0.5 


A = -0.6 


A = 1.5 


A = 2.0 


O.OO(exact) 


0.5424 


0.5470 


0.5513 


0.5553 


0.5589 


0.5619 


0.5641 


***** 


***** 


0.00 


0.5424 


0.5469 


0.5512 


0.5552 


0.5588 


0.5618 


0.5640 


***** 


***** 


0.01 


0.5423 


0.5469 


0.5512 


0.5551 


0.5587 


0.5617 


0.5640 


0.47(1) 


0.28(2) 


0.02 


0.5421 


0.5467 


0.5509 


0.5549 


0.5585 


0.5615 


0.5637 


0.453(4) 


0.350(6) 


0.03 


0.5418 


0.5463 


0.5506 


0.5545 


0.5581 


0.5611 


0.5633 


0.449(1) 


0.377(3) 


0.04 


0.5413 


0.5458 


0.5501 


0.5540 


0.5575 


0.5605 


0.5627 


0.448(1) 


0.390(2) 


0.05 


0.5407 


0.5452 


0.5494 


0.5533 


0.5568 


0.5598 


0.5620 


0.4482(9) 


0.399(1) 


0.06 


0.5400 


0.5444 


0.5486 


0.5525 


0.5560 


0.5589 


0.5610 


0.4493(6) 


0.406(1) 


0.07 


0.5391 


0.5435 


0.5477 


0.5515 


0.5549 


0.5578 


0.5599 


0.4505(6) 


0.4120(9) 


0.08 


0.5381 


0.5425 


0.5466 


0.5504 


0.5538 


0.5566 


0.5587 


0.4520(5) 


0.4172(7) 


0.09 


0.5370 


0.5413 


0.5453 


0.5491 


0.5524 


0.5552 


0.5572 


0.4536(4) 


0.4218(6) 


0.10 


0.5357 


0.5399 


0.5439 


0.5476 


0.5509 


0.5536 


0.5556 


0.4550(4) 


0.4258(5) 


0.11 


0.5343 


0.5384 


0.5423 


0.5460 


0.5492 


0.5519 


0.5538 


0.4563(3) 


0.4295(5) 


0.12 


0.5327 


0.5368 


0.5406 


0.5442 


0.5474 


0.5500 


0.5519 


0.4577(3) 


0.4329(4) 


0.13 


0.5310 


0.5350 


0.5388 


0.5423 


0.5454 


0.5479 


0.5497 


0.4589(3) 


0.4360(4) 


0.14 


0.5291 


0.5330 


0.5367 


0.5401 


0.5432 


0.5457 


0.5474 


0.4600(3) 


0.4387(4) 


0.15 


0.5271 


0.5309 


0.5345 


0.5379 


0.5408 


0.5432 


0.5449 


0.4609(3) 


0.4412(4) 


0.16 


0.5249 


0.5287 


0.5322 


0.5354 


0.5383 


0.5406 


0.5422 


0.4617(2) 


0.4435(3) 


0.17 


0.5226 


0.5263 


0.5297 


0.5328 


0.5356 


0.5378 


0.5393 


0.4623(3) 


0.4454(4) 


0.18 


0.5201 


0.5237 


0.5270 


0.5300 


0.5327 


0.5348 


0.5362 


0.4628(3) 


0.4471(3) 


0.19 


0.5175 


0.5209 


0.5241 


0.5270 


0.5296 


0.5316 


0.5329 


0.4631(2) 


0.4486(3) 


0.20 


0.5147 


0.5180 


0.5210 


0.5239 


0.5263 


0.5282 


0.5294 


0.4632(2) 


0.4497(3) 


0.21 


0.5117 


0.5148 


0.5178 


0.5205 


0.5228 


0.5247 


0.5257 


0.4630(2) 


0.4506(3) 


0.22 


0.5085 


0.5115 


0.5144 


0.5169 


0.5192 


0.5209 


0.5218 


0.4627(2) 


0.4513(3) 


0.23 


0.5052 


0.5080 


0.5107 


0.5132 


0.5153 


0.5169 


0.5177 


0.4622(2) 


0.4517(3) 


0.24 


0.5016 


0.5043 


0.5069 


0.5092 


0.5112 


0.5126 


0.5133 


0.4614(2) 


0.4517(3) 


0.25 


0.4979 


0.5004 


0.5029 


0.5051 


0.5069 


0.5082 


0.5087 


0.4604(2) 


0.4516(3) 


0.26 


0.4939 


0.4963 


0.4986 


0.5007 


0.5023 


0.5035 


0.5039 


0.4591(2) 


0.4511(2) 


0.27 


0.4897 


0.4920 


0.4941 


0.4960 


0.4976 


0.4986 


0.4988 


0.4575(2) 


0.4502(3) 


0.28 


0.4853 


0.4874 


0.4894 


0.4911 


0.4925 


0.4934 


0.4935 


0.4557(2) 


0.4491(3) 


0.29 


0.4806 


0.4826 


0.4844 


0.4860 


0.4873 


0.4880 


0.4879 


0.4536(2) 


0.4476(2) 


0.30 


0.4757 


0.4775 


0.4792 


0.4806 


0.4817 


0.4823 


0.4820 


0.4511(2) 


0.4457(2) 


0.31 


0.4705 


0.4721 


0.4736 


0.4749 


0.4759 


0.4763 


0.4758 


0.4483(2) 


0.4435(2) 


0.32 


0.4650 


0.4664 


0.4678 


0.4689 


0.4697 


0.4699 


0.4693 


0.4451(2) 


0.4409(2) 


0.33 


0.4591 


0.4605 


0.4617 


0.4626 


0.4632 


0.4633 


0.4625 


0.4415(2) 


0.4378(2) 


0.34 


0.4530 


0.4541 


0.4552 


0.4560 


0.4564 


0.4563(1) 


0.4553(1) 


0.4375(2) 


0.4343(2) 


0.35 


0.4464 


0.4474 


0.4483 


0.4489 


0.4492(1) 


0.4490(1) 


0.4478(1) 


0.4330(2) 


0.4303(2) 


0.36 


0.4394 


0.4403 


0.4410 


0.4415 


0.4416(1) 


0.4412(1) 


0.4399(1) 


0.4280(2) 


0.4257(2) 


0.37 


0.4320 


0.4327 


0.4332 


0.4336 


0.4336(1) 


0.4330(1) 


0.4315(1) 


0.4224(2) 


0.4205(2) 


0.38 


0.4241 


0.4246 


0.4250 


0.4252(1) 


0.4250(1) 


0.4243(1) 


0.4226(1) 


0.4163(2) 


0.4147(2) 


0.39 


0.4155 


0.4159 


0.4162 


0.4162(1) 


0.4159(1) 


0.4151(1) 


0.4132(1) 


0.4094(1) 


0.4081(2) 


0.40 


0.4063 


0.4066 


0.4067 


0.4066(1) 


0.4062(1) 


0.4052(1) 


0.4032(2) 


0.4017(1) 


0.4007(2) 


0.41 


0.3963 


0.3964 


0.3964 


0.3962(1) 


0.3957(1) 


0.3946(2) 


0.3925(2) 


0.3931(1) 


0.3923(2) 


0.42 


0.3854 


0.3854 


0.3853(1) 


0.3850(1) 


0.3844(1) 


0.3832(2) 


0.3810(2) 


0.3834(1) 


0.3829(1) 


0.43 


0.3734 


0.3733 


0.3731(1) 


0.3727(1) 


0.3720(2) 


0.3707(2) 


0.3685(3) 


0.3725(1) 


0.3722(1) 


0.44 


0.3600 


0.3598(1) 


0.3595(1) 


0.3590(1) 


0.3583(2) 


0.3570(2) 


0.3548(3) 


0.3599(1) 


0.3598(1) 


0.45 


0.3449 


0.3446(1) 


0.3442(1) 


0.3437(2) 


0.3429(2) 


0.3416(3) 


0.3395(4) 


0.3455(1) 


0.3455(1) 


0.46 


0.3274 


0.3271 


0.3267(1) 


0.3261(1) 


0.3253(2) 


0.3241(3) 


0.3220(4) 


0.3286(1) 


0.3286(1) 


0.47 


0.3062(2) 


0.3059(1) 


0.3055(1) 


0.3050 


0.3042 


0.3031(1) 


0.3013(2) 


0.3076(4) 


0.3077(4) 


0.48 


0.2791(2) 


0.2788(2) 


0.2785(2) 


0.2781(3) 


0.2774(3) 


0.2765(4) 


0.2751(6) 


0.2804(1) 


0.2806(1) 


0.49 


0.240(4) 


0.240(4) 


0.240(4) 


0.239(4) 


0.239(4) 


0.238(4) 


0.238(4) 


0.241(4) 


0.241(4) 
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(b) bi (1 > A > 0) 



M A = 1.0 


A = 0.9 


A = 0.8 


A = 0.7 


A = 0.6 


A = 0.5 


A = 0.4 


A = 0.3 


A = 0.2 


A = 0.1 


O.OO(exact) ***** 


0.3321 


0.3323 


0.3297 


0.3254 


0.3195 


0.3124 


0.3040 


0.2943 


0.2834 



0.00 
0.01 
0.02 
0.03 
0.04 
0.05 
0.06 
0.07 
0.08 
0.09 
0.10 
0.11 
0.12 
0.13 
0.14 
0.15 
0.16 
0.17 
0.18 
0.19 
0.20 
0.21 
0.22 
0.23 
0.24 
0.25 
0.26 
0.27 
0.28 
0.29 
0.30 
0.31 
0.32 
0.33 
0.34 
0.35 
0.36 
0.37 
0.38 
0.39 
0.40 
0.41 
0.42 
0.43 
0.44 
0.45 
0.46 
0.47 
0.48 
0.49 



0.348(7) 
0.355(3) 
0.3580(9) 
0.361(2) 
0.366(1) 
0.367(2) 
0.3697(3) 
0.3717(4) 
0.374(1) 
0.3767(8) 
0.379(1) 
0.3823(7) 
0.385(1) 
0.388(1) 
0.391(1) 
0.394(1) 
0.3971(5) 
0.4003(9) 
0.404(1) 
0.4078(7) 
0.4126(6) 
0.4180(8) 
0.423(1) 
0.429(1) 
0.436(1) 
0.442(1) 
0.449(1) 
0.457(1) 
0.4657(9) 
0.4759(6) 
0.4878(4) 
0.5013(7) 
0.515(1) 
0.532(2) 
0.551(1) 
0.569(1) 
0.591(2) 
0.618(2) 
0.650(1) 
0.6882(9) 
0.735(1) 
0.792(3) 
0.860(5) 
0.953(5) 
1.077(2) 
1.222(5) 
1.46(3) 
1.89(3) 
2.9(2) 



0.325(5) 
0.343(5) 
0.348(2) 
0.3516(9) 
0.355(1) 
0.3588(9) 
0.360(1) 
0.3629(1) 
0.3649(2) 
0.367(1) 
0.3698(6) 
0.3721(8) 
0.3755(5) 
0.378(1) 
0.381(1) 
0.384(1) 
0.387(1) 
0.3906(5) 
0.3939(8) 
0.397(1) 
0.4017(6) 
0.4066(5) 
0.4121(7) 
0.417(1) 
0.424(1) 
0.431(1) 
0.436(1) 
0.444(1) 
0.452(1) 
0.4608(9) 
0.4712(5) 
0.4832(4) 
0.4968(7) 
0.511(1) 
0.528(2) 
0.547(1) 
0.565(1) 
0.588(2) 
0.615(2) 
0.647(1) 
0.6856(9) 
0.732(1) 
0.790(3) 
0.858(5) 
0.951(5) 
1.076(2) 
1.221(5) 
1.46(3) 
1.89(3) 
2.9(2) 



0.326(2) 
0.338(3) 
0.342(1) 
0.3450(7) 
0.348(1) 
0.3513(8) 
0.353(1) 
0.3556 
0.3576(1) 
0.360(1) 
0.3625(5) 
0.3649(6) 
0.3681(4) 
0.371(1) 
0.374(1) 
0.3772(9) 
0.380(1) 
0.3836(4) 
0.3870(6) 
0.3906(8) 
0.3951(5) 
0.4001(4) 
0.4057(6) 
0.411(1) 
0.418(1) 
0.4245(9) 
0.4306(9) 
0.4382(9) 
0.446(1) 
0.4555(8) 
0.4660(5) 
0.4783(4) 
0.4920(6) 
0.506(1) 
0.523(2) 
0.543(1) 
0.561(1) 
0.584(1) 
0.612(2) 
0.644(1) 
0.6828(8) 
0.730(1) 
0.787(3) 
0.855(5) 
0.949(5) 
1.074(2) 
1.220(5) 
1.46(3) 
1.89(3) 
2.9(2) 



0.325(1) 
0.332(2) 
0.336(1) 
0.3382(6) 
0.341(1) 
0.3434(7) 
0.345(1) 
0.3479 
0.3498(1) 
0.3521(8) 
0.3546(4) 
0.3571(5) 
0.3603(3) 
0.3628(7) 
0.366(1) 
0.3694(8) 
0.372(1) 
0.3761(3) 
0.3796(5) 
0.3833(7) 
0.3879(4) 
0.3931(3) 
0.3987(5) 
0.404(1) 
0.411(1) 
0.4179(7) 
0.4242(8) 
0.4320(8) 
0.4405(9) 
0.4497(8) 
0.4604(4) 
0.4728(3) 
0.4867(6) 
0.501(1) 
0.518(1) 
0.538(1) 
0.557(1) 
0.580(1) 
0.608(2) 
0.640(1) 
0.6796(8) 
0.727(1) 
0.784(3) 
0.853(5) 
0.947(4) 
1.072(2) 
1.219(5) 
1.46(3) 
1.89(3) 
2.9(2) 



0.3222(4) 
0.326(1) 
0.3287(8) 
0.3309(5) 
0.3329(9) 
0.3352(6) 
0.3374(7) 
0.3396(1) 

0.3415 
0.3437(7) 
0.3462(3) 
0.3487(3) 
0.3518(2) 
0.3544(6) 
0.358(1) 
0.3610(6) 
0.3642(8) 
0.3680(2) 
0.3715(4) 
0.3755(6) 
0.3801(3) 
0.3854(2) 
0.3912(4) 
0.3967(9) 
0.403(1) 
0.4107(6) 
0.4173(7) 
0.4252(7) 
0.4339(9) 
0.4433(7) 
0.4543(4) 
0.4669(3) 
0.4809(5) 
0.495(1) 
0.513(1) 
0.5327(9) 
0.5519(9) 
0.576(1) 
0.604(2) 
0.636(1) 
0.6762(8) 
0.723(1) 
0.781(3) 
0.850(5) 
0.944(4) 
1.070(2) 
1.218(5) 
1.46(3) 
1.89(3) 
2.9(2) 



0.3173(1) 
0.320(1) 
0.3214(6) 
0.3232(3) 
0.3247(8) 
0.3266(5) 
0.3287(5) 
0.3309(1) 
0.3326(1) 
0.3348(5) 
0.3372(2) 
0.3397(2) 
0.3428(1) 
0.3454(4) 
0.3485(8) 
0.3520(5) 
0.3553(6) 
0.3592(1) 
0.3628(3) 
0.3670(5) 
0.3717(2) 
0.3771(1) 
0.3829(3) 
0.3886(7) 
0.396(1) 
0.4029(5) 
0.4097(6) 
0.4178(6) 
0.4267(8) 
0.4364(6) 
0.4476(3) 
0.4604(3) 
0.4745(5) 
0.489(1) 
0.507(1) 
0.5269(8) 
0.5466(8) 
0.571(1) 
0.599(2) 
0.632(1) 
0.6724(8) 
0.720(1) 
0.778(3) 
0.847(5) 
0.941(4) 
1.067(2) 
1.217(6) 
1.46(3) 
1.89(3) 
2.9(2) 



0.3109 
0.3121(7) 
0.3133(4) 
0.3147(2) 
0.3159(6) 
0.3174(4) 
0.3195(4) 
0.3215(1) 
0.3232(1) 
0.3253(4) 
0.3276(2) 
0.3301(1) 

0.3330 
0.3357(3) 
0.3388(6) 
0.3423(4) 
0.3457(4) 
0.3496(1) 
0.3535(2) 
0.3577(4) 
0.3626(2) 
0.3681(1) 
0.3740(2) 
0.3798(6) 
0.3868(8) 
0.3943(4) 
0.4014(4) 
0.4098(5) 
0.4188(7) 
0.4288(6) 
0.4403(3) 
0.4532(2) 
0.4674(5) 
0.482(1) 
0.500(1) 
0.5205(7) 
0.5407(6) 
0.565(1) 
0.594(2) 
0.628(1) 
0.6682(8) 
0.716(1) 
0.774(3) 
0.843(5) 
0.938(3) 
1.065(2) 
1.215(6) 
1.46(3) 
1.89(3) 
2.9(2) 



0.3031 
0.3036(5) 
0.3044(3) 
0.3055(1) 
0.3065(5) 
0.3077(4) 
0.3096(3) 
0.3114(1) 
0.3130(1) 
0.3151(3) 
0.3173(1) 
0.3198(1) 

0.3226 
0.3253(2) 
0.3284(5) 
0.3319(3) 
0.3354(2) 

0.3394 
0.3433(2) 
0.3477(4) 
0.3527(1) 

0.3582 
0.3642(1) 
0.3702(4) 
0.3773(6) 
0.3849(3) 
0.3923(3) 
0.4009(4) 
0.4102(6) 
0.4205(5) 
0.4321(2) 
0.4452(1) 
0.4596(4) 
0.4748(9) 
0.493(1) 
0.5133(5) 
0.5342(5) 
0.559(1) 
0.588(1) 
0.622(1) 
0.664(1) 
0.711(1) 
0.769(4) 
0.839(4) 
0.934(3) 
1.061(2) 
1.214(6) 

1.46(3) 

1.89(3) 

2.9(3) 



0.2938 
0.2941(3) 
0.2946(2) 

0.2954 
0.2962(4) 
0.2972(3) 
0.2989(2) 
0.3006(1) 
0.3022(1) 
0.3041(2) 
0.3062(1) 
0.3087 
0.3114 
0.3141(1) 
0.3172(3) 
0.3206(2) 
0.3243(1) 

0.3282 
0.3323(1) 
0.3368(3) 
0.3419(1) 

0.3475 
0.3535(1) 
0.3597(3) 
0.3668(5) 
0.3745(2) 
0.3823(2) 
0.3911(3) 
0.4006(5) 
0.4112(4) 
0.4232(1) 
0.4363(1) 
0.4508(4) 
0.4664(7) 
0.4847(9) 
0.5053(4) 
0.5268(3) 
0.552(1) 
0.581(1) 
0.617(1) 
0.658(1) 
0.706(1) 
0.764(4) 
0.835(5) 
0.930(3) 
1.058(2) 
1.212(6) 
1.46(3) 
1.89(3) 
3.0(3) 



0.2832 
0.2833(1) 
0.2837(1) 

0.2843 
0.2850(2) 
0.2859(2) 
0.2873(1) 
0.2889 
0.2904 
0.2923(1) 
0.2943(1) 
0.2967 
0.2993 
0.3020 
0.3051(2) 
0.3085(1) 
0.3122 
0.3162 
0.3203 
0.3250(2) 
0.3301 
0.3357 
0.3418(1) 
0.3481(2) 
0.3553(3) 
0.3631(1) 
0.3713(1) 
0.3803(2) 
0.3901(3) 
0.4010(3) 
0.4132(1) 
0.4264(1) 
0.4410(4) 
0.4569(6) 
0.4754(7) 
0.4962(3) 
0.5186(1) 
0.544(1) 
0.574(1) 
0.610(1) 
0.652(1) 
0.700(2) 
0.759(4) 
0.830(5) 
0.925(2) 
1.054(2) 
1.210(6) 
1.46(3) 
1.89(3) 
3.0(3) 
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(b) bi (0 > A > -0.6, A = 1.5, 2.0) 



M 


A = 0.0 


A = -0.1 


A = -0.2 


A = -0.3 


A = -0.4 


A = -0.5 


A = -0.6 


A = 1.5 


A = 2.0 


O.OO(exact) 


0.2712 


0.2576 


0.2424 


0.2255 


0.2068 


0.1861 


0.1630 


***** 


***** 


0.00 


0.2712 


0.2576 


0.2425 


0.2255(1) 


0.2063(4) 


0.1844(9) 


0.158(2) 


***** 


***** 


0.01 


0.2712 


0.2577(1) 


0.2425(3) 


0.2255(7) 


0.206(1) 


0.184(4) 


0.16(1) 


0.37(2) 


0.279(4) 


0.02 


0.2716 


0.2581(1) 


0.2432(3) 


0.2266(7) 


0.208(1) 


0.188(2) 


0.166(7) 


0.39(1) 


0.38(1) 


0.03 


0.2721 


0.2585(1) 


0.2435(3) 


0.2269(7) 


0.209(1) 


0.188(3) 


0.166(7) 


0.3867(1) 


0.395(3) 


0.04 


0.2728 


0.2593(2) 


0.2446(6) 


0.228(1) 


0.211(2) 


0.191(4) 


0.171(8) 


0.388(6) 


0.40(1) 


0.05 


0.2737 


0.2603(1) 


0.2458(5) 


0.230(1) 


0.213(1) 


0.195(3) 


0.179(6) 


0.396(2) 


0.411(4) 


0.06 


0.2748 


0.2611(1) 


0.2461(4) 


0.230(1) 


0.212(2) 


0.192(4) 


0.170(9) 


0.394(5) 


0.408(7) 


0.07 


0.2761 


0.2622 


0.2469(1) 


0.2300(1) 


0.2112(2) 


0.1898(4) 


0.165(1) 


0.397(1) 


0.412(3) 


0.08 


0.2777 


0.2638(1) 


0.2486(4) 


0.2317(9) 


0.213(2) 


0.192(4) 


0.17(1) 


0.399(1) 


0.415(3) 


0.09 


0.2794 


0.2654 


0.2500(1) 


0.2330(4) 


0.214(1) 


0.192(2) 


0.165(7) 


0.401(3) 


0.417(5) 


0.10 


0.2814 


0.2674(1) 


0.2520(1) 


0.2349(3) 


0.2157(8) 


0.193(1) 


0.165(4) 


0.404(1) 


0.421(2) 


0.11 


0.2837 


0.2697 


0.2544(1) 


0.2376(3) 


0.2192(8) 


0.199(1) 


0.177(4) 


0.406(2) 


0.423(3) 


0.12 


0.2862 


0.2721 


0.2568 


0.2402(1) 


0.2221(3) 


0.2026(9) 


0.182(2) 


0.409(1) 


0.427(2) 


0.13 


0.2890 


0.2749 


0.2596(2) 


0.2429(7) 


0.225(1) 


0.205(4) 


0.18(1) 


0.412(2) 


0.429(3) 


0.14 


0.2920(1) 


0.2779 


0.2626(3) 


0.2461(8) 


0.228(1) 


0.209(3) 


0.190(8) 


0.414(3) 


0.432(4) 


0.15 


0.2954 


0.2812 


0.2659(1) 


0.2493(4) 


0.2314(9) 


0.212(2) 


0.193(4) 


0.418(2) 


0.436(3) 


0.16 


0.2991 


0.2848 


0.2693 


0.2523 


0.2335(2) 


0.2125(8) 


0.188(2) 


0.419(3) 


0.437(4) 


0.17 


0.3031 


0.2888 


0.2732(1) 


0.2561(3) 


0.2372(6) 


0.216(1) 


0.191(3) 


0.423(1) 


0.441(1) 


0.18 


0.3074 


0.2933 


0.2778(1) 


0.2609(4) 


0.242(1) 


0.221(2) 


0.197(7) 


0.426(1) 


0.443(2) 


0.19 


0.3121(1) 


0.2981(1) 


0.2827(2) 


0.2658(4) 


0.2471(8) 


0.226(1) 


0.201(5) 


0.428(1) 


0.446(2) 


0.20 


0.3173(1) 


0.3033(2) 


0.2880(3) 


0.2712(6) 


0.253(1) 


0.232(1) 


0.208(4) 


0.432(1) 


0.449(1) 


0.21 


0.3229 


0.3089(1) 


0.2934(2) 


0.2765(4) 


0.2579(7) 


0.238(1) 


0.216(1) 


0.437(1) 


0.453(1) 


0.22 


0.3289(1) 


0.3148(1) 


0.2993(2) 


0.2823(2) 


0.2635(1) 


0.2430(3) 


0.221(1) 


0.442(1) 


0.458(1) 


0.23 


0.3354(1) 


0.3215(3) 


0.3061(4) 


0.2891(6) 


0.2704(8) 


0.250(1) 


0.227(3) 


0.446(2) 


0.463(2) 


0.24 


0.3426(2) 


0.3286(2) 


0.3131(3) 


0.2959(6) 


0.277(1) 


0.256(2) 


0.232(5) 


0.452(2) 


0.468(3) 


0.25 


0.3505(1) 


0.3365(1) 


0.3210(2) 


0.3038(4) 


0.2845(8) 


0.263(1) 


0.238(4) 


0.458(1) 


0.474(2) 


0.26 


0.3590 


0.3455(1) 


0.3304(3) 


0.3136(6) 


0.2949(9) 


0.274(1) 


0.250(1) 


0.463(1) 


0.478(2) 


0.27 


0.3683(1) 


0.3550 


0.3402(1) 


0.3238(2) 


0.3054(3) 


0.2849(5) 


0.2620(6) 


0.470(1) 


0.485(2) 


0.28 


0.3784(2) 


0.3654(2) 


0.3509(6) 


0.335(1) 


0.317(1) 


0.297(2) 


0.274(4) 


0.477(1) 


0.492(1) 


0.29 


0.3897(3) 


0.3771(5) 


0.3630(8) 


0.347(1) 


0.330(2) 


0.311(3) 


0.291(6) 


0.485(1) 


0.499(1) 


0.30 


0.4021(3) 


0.3897(5) 


0.3758(8) 


0.360(1) 


0.343(2) 


0.324(3) 


0.304(7) 


0.4947(7) 


0.5081(9) 


0.31 


0.4153(3) 


0.4028(5) 


0.3887(8) 


0.373(1) 


0.354(2) 


0.333(3) 


0.309(4) 


0.5060(5) 


0.5188(7) 


0.32 


0.4300(5) 


0.4175(7) 


0.403(1) 


0.387(1) 


0.368(2) 


0.346(3) 


0.321(5) 


0.519(1) 


0.531(1) 


0.33 


0.4462(7) 


0.434(1) 


0.420(1) 


0.404(2) 


0.386(4) 


0.364(6) 


0.34(1) 


0.532(2) 


0.544(2) 


0.34 


0.4648(5) 


0.4527(7) 


0.439(1) 


0.422(2) 


0.403(3) 


0.380(5) 


0.351(9) 


0.548(2) 


0.560(3) 


0.35 


0.4859(2) 


0.4741(3) 


0.4605(6) 


0.445(1) 


0.426(2) 


0.404(3) 


0.376(7) 


0.567(1) 


0.577(1) 


0.36 


0.5092(1) 


0.4986(4) 


0.4864(8) 


0.472(1) 


0.457(2) 


0.438(3) 


0.418(4) 


0.583(1) 


0.593(1) 


0.37 


0.5356(8) 


0.5257(5) 


0.5144(4) 


0.5015(7) 


0.487(1) 


0.470(3) 


0.453(6) 


0.605(2) 


0.614(2) 


0.38 


0.566(1) 


0.557(1) 


0.546(1) 


0.534(2) 


0.520(4) 


0.504(7) 


0.49(1) 


0.631(2) 


0.640(2) 


0.39 


0.603(1) 


0.594(2) 


0.585(3) 


0.574(4) 


0.563(6) 


0.550(9) 


0.54(1) 


0.661(1) 


0.669(1) 


0.40 


0.646(2) 


0.638(2) 


0.629(3) 


0.620(4) 


0.608(6) 


0.60(1) 


0.58(1) 


0.699(1) 


0.706(1) 


0.41 


0.694(2) 


0.686(3) 


0.677(4) 


0.666(5) 


0.653(7) 


0.64(1) 


0.62(1) 


0.744(1) 


0.751(1) 


0.42 


0.752(5) 


0.744(5) 


0.735(6) 


0.724(8) 


0.71(1) 


0.69(1) 


0.67(2) 


0.801(3) 


0.808(3) 


0.43 


0.823(5) 


0.816(6) 


0.808(8) 


0.80(1) 


0.78(1) 


0.77(2) 


0.75(3) 


0.868(6) 


0.874(6) 


0.44 


0.919(3) 


0.912(4) 


0.904(5) 


0.893(7) 


0.88(1) 


0.86(1) 


0.84(2) 


0.961(6) 


0.966(7) 


0.45 


1.049(2) 


1.043(3) 


1.036(3) 


1.028(5) 


1.018(7) 


1.01(1) 


0.99(1) 


1.084(2) 


1.088(2) 


0.46 


1.208(6) 


1.206(7) 


1.203(7) 


1.200(7) 


1.198(8) 


1.196(9) 


1.20(1) 


1.225(5) 


1.228(5) 


0.47 


1.46(3) 


1.46(3) 


1.46(3) 


1.46(3) 


1.46(3) 


1.46(3) 


1.46(3) 


1.46(3) 


1.47(3) 


0.48 


1.89(3) 


1.89(3) 


1.89(3) 


1.89(2) 


1.89(2) 


1.89(2) 


1.89(3) 


1.89(4) 


1.89(4) 


0.49 


3.0(3) 


3.0(3) 


3.0(3) 


3.0(3) 


3.0(3) 


3.0(3) 


3.0(3) 


2.9(2) 


2.9(2) 
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(c) ai (1 > A > 0) 



M A = 1.0 


A = 0.9 


A = 0.8 


A = 0.7 


A = 0.6 


A = 0.5 


A = 0.4 


A = 0.3 


A = 0.2 


A = 0.1 


O.OO(exact) ***** 


0.7049 


0.6069 


0.5464 


0.5008 


0.4629 


0.4297 


0.3995 


0.3712 


0.3444 



0.00 
0.01 
0.02 
0.03 
0.04 
0.05 
0.06 
0.07 
0.08 
0.09 
0.10 
0.11 
0.12 
0.13 
0.14 
0.15 
0.16 
0.17 
0.18 
0.19 
0.20 
0.21 
0.22 
0.23 
0.24 
0.25 
0.26 
0.27 
0.28 
0.29 
0.30 
0.31 
0.32 
0.33 
0.34 
0.35 
0.36 
0.37 
0.38 
0.39 
0.40 
0.41 
0.42 
0.43 
0.44 
0.45 
0.46 
0.47 
0.48 
0.49 



0.626(8) 
0.590(2) 
0.568(1) 
0.549(5) 
0.5375(6) 
0.520(3) 
0.506(1) 
0.4947(6) 
0.484(3) 
0.475(1) 
0.464(2) 
0.454(1) 
0.4456(7) 
0.437(2) 
0.4298(9) 
0.422(1) 
0.4140(9) 
0.4070(7) 
0.400(1) 
0.3938(7) 
0.387(1) 
0.3813(6) 
0.3756(5) 
0.3700(7) 
0.3644(5) 
0.3596(7) 
0.3545(2) 
0.3499(4) 
0.3455(4) 
0.3406(3) 
0.3371(4) 
0.3330(2) 
0.3294(3) 
0.3262(1) 
0.3221(2) 
0.3200(3) 
0.3169(5) 
0.3145(4) 

0.3124 
0.3094(1) 
0.3088(3) 
0.3070(8) 
0.3062(6) 
0.3058(1) 
0.3045(1) 
0.3062(5) 
0.307(1) 
0.3094(9) 
0.313(1) 



0.64(1) 
0.589(6) 
0.562(2) 
0.545(1) 
0.529(4) 
0.5190(4) 
0.504(3) 
0.4921(9) 
0.4818(4) 
0.472(2) 
0.4641(8) 
0.454(2) 
0.4453(9) 
0.4375(6) 
0.430(1) 
0.4229(7) 
0.415(1) 
0.4083(8) 
0.4019(5) 
0.396(1) 
0.3895(6) 
0.384(1) 
0.3778(5) 
0.3725(4) 
0.3673(7) 
0.3619(4) 
0.3575(6) 
0.3526(2) 
0.3483(3) 
0.3441(3) 
0.3394(3) 
0.3362(4) 
0.3323(2) 
0.3289(3) 
0.3258(1) 
0.3219(2) 
0.3199(3) 
0.3170(4) 
0.3147(4) 

0.3128 
0.3098(1) 
0.3093(3) 
0.3076(7) 
0.3068(6) 
0.3065(1) 
0.3052(1) 
0.3069(5) 
0.308(1) 
0.3099(9) 
0.313(1) 



0.587(6) 
0.553(4) 
0.534(1) 
0.520(1) 
0.508(3) 
0.4995(2) 
0.487(2) 
0.4767(7) 
0.4678(3) 
0.459(2) 
0.4522(5) 
0.443(1) 
0.4356(7) 
0.4285(4) 
0.421(1) 
0.4153(6) 
0.408(1) 
0.4021(6) 
0.3962(4) 
0.390(1) 
0.3847(5) 
0.3793(8) 
0.3740(4) 
0.3690(3) 
0.3642(6) 
0.3592(4) 
0.3550(5) 
0.3505(1) 
0.3465(3) 
0.3426(3) 
0.3382(2) 
0.3352(3) 
0.3315(2) 
0.3283(3) 
0.3254(1) 
0.3217(1) 
0.3199(3) 
0.3171(4) 
0.3149(3) 

0.3131 
0.3103(1) 
0.3098(3) 
0.3083(7) 
0.3075(5) 
0.3072(1) 
0.3060(1) 
0.3076(5) 
0.308(1) 
0.3104(8) 
0.3135(9) 



0.540(2) 
0.519(3) 
0.505(1) 
0.4953(9) 
0.486(2) 
0.4790(1) 
0.468(2) 
0.4602(5) 
0.4527(2) 
0.445(1) 
0.4394(4) 
0.432(1) 
0.4249(6) 
0.4187(3) 
0.412(1) 
0.4069(4) 
0.401(1) 
0.3951(5) 
0.3898(3) 
0.3845(8) 
0.3794(4) 
0.3745(7) 
0.3696(3) 
0.3651(3) 
0.3607(5) 
0.3561(3) 
0.3523(4) 
0.3482(1) 
0.3444(2) 
0.3408(2) 
0.3367(2) 
0.3340(3) 
0.3306(2) 
0.3276(2) 
0.3250(1) 
0.3215(1) 
0.3198(2) 
0.3172(4) 
0.3152(3) 

0.3135 
0.3109(1) 
0.3105(3) 
0.3090(7) 
0.3083(5) 
0.3080(1) 
0.3068(1) 
0.3083(5) 
0.3089(9) 
0.3110(8) 
0.3138(9) 



0.4990(5) 
0.486(2) 
0.4771(7) 
0.4699(7) 
0.463(2) 

0.4577 
0.449(1) 
0.4427(4) 
0.4366(2) 
0.430(1) 
0.4255(2) 
0.419(1) 
0.4132(4) 
0.4078(2) 
0.402(1) 
0.3977(3) 
0.3923(8) 
0.3873(3) 
0.3827(2) 
0.3780(7) 
0.3735(3) 
0.3692(6) 
0.3648(2) 
0.3607(2) 
0.3568(4) 
0.3526(2) 
0.3492(4) 

0.3455 
0.3421(2) 
0.3389(2) 
0.3351(1) 
0.3327(3) 
0.3296(2) 
0.3269(2) 

0.3244 
0.3212(1) 
0.3197(2) 
0.3174(4) 
0.3155(3) 

0.3140 
0.3116(1) 
0.3112(3) 
0.3098(6) 
0.3092(5) 
0.3089(1) 
0.3078(1) 
0.3092(4) 
0.3097(8) 
0.3116(7) 
0.3142(8) 



0.4626(1) 
0.455(1) 
0.4493(5) 
0.4444(5) 
0.439(1) 

0.4356 
0.429(1) 
0.4241(3) 
0.4193(1) 
0.414(1) 
0.4104(1) 
0.4051(9) 
0.4004(3) 
0.3960(1) 
0.3915(8) 
0.3874(2) 
0.3829(7) 
0.3787(2) 
0.3747(1) 
0.3707(5) 
0.3668(2) 
0.3631(5) 
0.3593(1) 
0.3558(1) 
0.3523(3) 
0.3487(1) 
0.3457(3) 

0.3424 
0.3394(1) 
0.3366(1) 
0.3333(1) 
0.3311(2) 
0.3283(2) 
0.3260(2) 

0.3238 
0.3209(1) 
0.3196(2) 
0.3175(4) 
0.3158(3) 

0.3145 

0.3123 
0.3120(2) 
0.3108(6) 
0.3102(4) 
0.3099(1) 
0.3090(1) 
0.3102(4) 
0.3107(7) 
0.3123(6) 
0.3147(7) 



0.4297(4) 
0.425(1) 
0.4220(3 
0.4189(3 
0.415(1) 

0.4129 
0.4083(9 
0.4046(2 
0.4010(1 
0.3973(8 

0.3943 
0.3901(7 
0.3864(2 

0.3830 
0.3794(6 
0.3761(1 
0.3725(5 
0.3691(1 

0.3658 
0.3625(4 
0.3593(1 
0.3563(3 

0.3531 
0.3502(1 
0.3473(2 
0.3442(1 
0.3417(2 

0.3389 
0.3364(1 
0.3340(1 

0.3311 
0.3293(2 
0.3269(2 
0.3249(2 

0.3231 

0.3206 
0.3195(2 
0.3176(3 
0.3162(3 

0.3150 

0.3131 
0.3129(2 
0.3118(5 
0.3113(4 
0.3111(1 

0.3103 
0.3114(3 
0.3117(6 
0.3132(5 
0.3152(6 



0.3995(3) 
0.3971(8) 
0.3953(2) 
0.3935(1) 
0.3913(8) 

0.3897 
0.3867(6) 
0.3842(1) 
0.3817(1) 
0.3791(6) 

0.3769 
0.3740(4) 
0.3713(1) 

0.3687 
0.3661(4) 

0.3637 
0.3610(3) 

0.3583 

0.3559 
0.3534(3) 

0.3508 
0.3485(2) 

0.3460 

0.3438 
0.3415(2) 

0.3390 
0.3371(2) 
0.3349(1) 
0.3329(1) 
0.3310(1) 

0.3286 
0.3272(1) 
0.3253(2) 
0.3237(1) 

0.3222 

0.3201 
0.3193(1) 
0.3178(3) 
0.3166(2) 

0.3157 

0.3141 
0.3139(2) 
0.3131(4) 
0.3127(3) 
0.3125(1) 

0.3118 
0.3127(3) 
0.3130(5) 
0.3142(4) 
0.3158(5) 



0.3713(2 
0.3701(4 
0.3692(1 
0.3683(1 
0.3670(5 
0.3662 
0.3644(4 
0.3629(1 
0.3615 
0.3599(3 
0.3585 
0.3566(3 
0.3549 
0.3533 
0.3515(2 
0.3499 
0.3482(2 
0.3464 
0.3447 
0.3430(2 
0.3413 
0.3397(1 
0.3380 
0.3364 
0.3349(1 
0.3331 
0.3318(1 
0.3302(1 
0.3288(1 
0.3274 
0.3257 
0.3248(1 
0.3234(1 
0.3222(1 
0.3211 
0.3196 
0.3190(1 
0.3179(2 
0.3171(1 
0.3164 
0.3153 
0.3151(1 
0.3145(3 
0.3142(2 
0.3141 
0.3136 
0.3142(2 
0.3145(4 
0.3153(3 
0.3165(3 
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(c) ai (0 > A > -0.6, A = 1.5, 2.0) 



M 


A = 0.0 


A = -0.1 


A = -0.2 


A = -0.3 


A = -0.4 


A = -0.5 


A = -0.6 


A = 1.5 


A = 2.0 


O.OO(exact) 


0.3183 


0.2928 


0.2674 


0.2418 


0.2160 


0.1895 


0.1621 


***** 


***** 


0.00 


0.3183 


0.2927 


0.2674(2) 


0.2420(2) 


0.2162(4) 


0.1895(7) 


0.161(2) 


***** 


***** 


0.01 


0.3183 


0.2930(1) 


0.2678(2) 


0.2425(3) 


0.2167(4) 


0.1903(4) 


0.1630(5) 


0.80(2) 


0.92(5) 


0.02 


0.3183 


0.2932 


0.2681(1) 


0.2428(1) 


0.2171(1) 


0.1907(2) 


0.1633(2) 


0.712(6) 


0.789(8) 


0.03 


0.3183 


0.2934 


0.2684 


0.2432 


0.2174 


0.1910 


0.1637 


0.666(2) 


0.728(3) 


0.04 


0.3183 


0.2937(2) 


0.2690(3) 


0.2438(5) 


0.2181(7) 


0.1917(9) 


0.164(1) 


0.63(1) 


0.69(1) 


0.05 


0.3183 


0.2940 


0.2693 


0.2442 


0.2186 


0.1925 


0.1658(2) 


0.612(2) 


0.660(5) 


0.06 


0.3183 


0.2945(1) 


0.2702(3) 


0.2454(4) 


0.2198(6) 


0.1933(8) 


0.165(1) 


0.583(7) 


0.62(1) 


0.07 


0.3183 


0.2950 


0.2711 


0.2464 


0.2211(1) 


0.1948(3) 


0.168(1) 


0.562(2) 


0.598(4) 


0.08 


0.3183 


0.2955 


0.2719 


0.2475(1) 


0.2223(2) 


0.1961(4) 


0.169(1) 


0.545(1) 


0.577(3) 


0.09 


0.3183 


0.2961(1) 


0.2730(3) 


0.2489(5) 


0.2237(8) 


0.197(1) 


0.169(2) 


0.529(5) 


0.558(7) 


0.10 


0.3183 


0.2966 


0.2739 


0.2503 


0.2257(1) 


0.2005(3) 


0.1757(9) 


0.517(2) 


0.543(3) 


0.11 


0.3183 


0.2973(1) 


0.2752(3) 


0.2518(5) 


0.2270(8) 


0.201(1) 


0.172(3) 


0.501(4) 


0.524(5) 


0.12 


0.3183 


0.2981 


0.2765 


0.2536(2) 


0.2294(4) 


0.204(1) 


0.177(3) 


0.488(2) 


0.509(3) 


0.13 


0.3183 


0.2987 


0.2778(1) 


0.2553(3) 


0.2314(5) 


0.206(1) 


0.179(2) 


0.476(1) 


0.496(2) 


0.14 


0.3183 


0.2995(1) 


0.2791(3) 


0.2572(5) 


0.2335(9) 


0.208(1) 


0.179(4) 


0.465(3) 


0.483(4) 


0.15 


0.3183 


0.3003 


0.2807 


0.2597(1) 


0.2372(3) 


0.2139(7) 


0.192(1) 


0.456(1) 


0.472(2) 


0.16 


0.3183 


0.3010(1) 


0.2821(2) 


0.2612(5) 


0.238(1) 


0.213(2) 


0.184(6) 


0.445(2) 


0.459(3) 


0.17 


0.3183 


0.3019 


0.2838(1) 


0.2638(4) 


0.2419(9) 


0.218(2) 


0.192(6) 


0.435(1) 


0.448(2) 


0.18 


0.3183 


0.3027 


0.2854(2) 


0.2660(4) 


0.2446(9) 


0.221(2) 


0.195(5) 


0.426(1) 


0.438(1) 


0.19 


0.3183 


0.3036(1) 


0.2869(2) 


0.2683(5) 


0.247(1) 


0.223(2) 


0.195(6) 


0.418(2) 


0.429(2) 


0.20 


0.3183 


0.3046 


0.2891(1) 


0.2718(2) 


0.2527(4) 


0.233(1) 


0.215(2) 


0.410(1) 


0.420(1) 


0.21 


0.3183 


0.3053(1) 


0.2904(2) 


0.2733(6) 


0.254(1) 


0.231(3) 


0.202(9) 


0.402(1) 


0.411(2) 


0.22 


0.3183 


0.3062(1) 


0.2924(3) 


0.2765(7) 


0.258(1) 


0.238(3) 


0.215(9) 


0.394(1) 


0.402(1) 


0.23 


0.3183 


0.3071(1) 


0.2941(3) 


0.2790(6) 


0.262(1) 


0.241(2) 


0.217(7) 


0.3871(9) 


0.394(1) 


0.24 


0.3183 


0.3080 


0.2958(2) 


0.2815(5) 


0.265(1) 


0.244(2) 


0.218(7) 


0.380(1) 


0.387(1) 


0.25 


0.3183 


0.3091 


0.2983(1) 


0.2858(2) 


0.2716(5) 


0.256(1) 


0.243(2) 


0.3734(9) 


0.379(1) 


0.26 


0.3183 


0.3097(1) 


0.2995(3) 


0.2872(7) 


0.272(1) 


0.254(3) 


0.23(1) 


0.368(1) 


0.372(1) 


0.27 


0.3183 


0.3107(1) 


0.3017(4) 


0.2908(9) 


0.278(2) 


0.262(4) 


0.24(1) 


0.3614(5) 


0.3657(8) 


0.28 


0.3183 


0.3115(1) 


0.3034(3) 


0.2935(8) 


0.281(1) 


0.266(3) 


0.247(9) 


0.3558(6) 


0.3595(8) 


0.29 


0.3183 


0.3123 


0.3051(1) 


0.2961(4) 


0.285(1) 


0.270(2) 


0.250(7) 


0.3505(6) 


0.3536(7) 


0.30 


0.3183 


0.3135 


0.3077 


0.3008(2) 


0.2927(4) 


0.2839(9) 


0.277(2) 


0.3447(5) 


0.3472(7) 


0.31 


0.3183 


0.3140(1) 


0.3087(3) 


0.3020(7) 


0.293(1) 


0.281(4) 


0.26(1) 


0.3405(5) 


0.3425(7) 


0.32 


0.3183 


0.3149(2) 


0.3108(5) 


0.306(1) 


0.299(2) 


0.291(5) 


0.28(1) 


0.3356(2) 


0.3372(2) 


0.33 


0.3183 


0.3157(1) 


0.3124(4) 


0.3082(9) 


0.303(1) 


0.295(4) 


0.285(9) 


0.3313(4) 


0.3324(5) 


0.34 


0.3183 


0.3164 


0.3139(1) 


0.3106(3) 


0.3062(9) 


0.300(2) 


0.290(6) 


0.3274(2) 


0.3282(3) 


0.35 


0.3183 


0.3174 


0.3164 


0.3152(1) 


0.3139(2) 


0.3133(6) 


0.316(1) 


0.3226(3) 


0.3229(4) 


0.36 


0.3183 


0.3178(1) 


0.3170(3) 


0.3159(7) 


0.314(1) 


0.311(3) 


0.306(9) 


0.3200(4) 


0.3200(4) 


0.37 


0.3183 


0.3185(2) 


0.3188(6) 


0.319(1) 


0.319(2) 


0.320(5) 


0.32(1) 


0.3165(5) 


0.3161(5) 


0.38 


0.3183 


0.3191(1) 


0.3199(4) 


0.3210(9) 


0.322(1) 


0.324(3) 


0.326(8) 


0.3136(4) 


0.3129(5) 


0.39 


0.3183 


0.3195 


0.3210(1) 


0.3228(3) 


0.3251(8) 


0.328(1) 


0.331(4) 


0.3112 


0.3103 


0.40 


0.3183 


0.3204 


0.3230 


0.3265 


0.3314 


0.3388(1) 


0.3516(4) 


0.3076(2) 


0.3064(2) 


0.41 


0.3183 


0.3204(1) 


0.3230(3) 


0.3263(7) 


0.331(1) 


0.336(2) 


0.344(6) 


0.3068(4) 


0.3056(4) 


0.42 


0.3183 


0.3209(2) 


0.3241(6) 


0.328(1) 


0.334(2) 


0.342(4) 


0.355(8) 


0.3048(9) 


0.3034(9) 


0.43 


0.3183 


0.3211(2) 


0.3245(5) 


0.329(1) 


0.335(1) 


0.343(3) 


0.356(6) 


0.3038(7) 


0.3023(7) 


0.44 


0.3183 


0.3211 


0.3247(2) 


0.3294(4) 


0.3356(7) 


0.344(1) 


0.358(3) 


0.3034(1) 


0.3018(1) 


a a c 
0.40 


a o i oo 
0.3183 


a 001 c 
U.3215 


0.32OO 


0.3311(1 J 


0.3350(2 ) 


A O A AK ( OA 

0.3490(3) 


A OC71 / A \ 

0.36/ 1(4) 


0.3018(2) 


0.3001(2) 


0.46 


0.3183 


0.3211(1) 


0.3245(4) 


0.3290(7) 


0.335(1) 


0.344(2) 


0.357(4) 


0.3039(6) 


0.3023(7) 


0.47 


0.3183 


0.3209(2) 


0.3242(6) 


0.329(1) 


0.334(2) 


0.343(3) 


0.356(5) 


0.305(1) 


0.303(1) 


0.48 


0.3183 


0.3203(2) 


0.3229(5) 


0.326(1) 


0.331(1) 


0.337(2) 


0.347(4) 


0.308(1) 


0.307(1) 


0.49 


0.3183 


0.3195(2) 


0.3211(5) 


0.323(1) 


0.326(1) 


0.330(2) 


0.335(4) 


0.312(1) 


0.311(1) 



